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Abstract 



^ , A qualitative mechanism for the emergence of domain structured background 

(~| , gluon fields due to singularities in gauge field configurations is considered, 

- , , and a model displaying a type of mean field approximation to the QCD par- 

X, 

. tition function based on this mechanism is formulated. Estimation of the 

vacuum parameters (gluon condensate, topological susceptibility, string con- 
stant and quark condensate) indicates that domain-like structures lead to an 
area law for the Wilson loop, nonzero topological susceptibility and sponta- 
neous breakdown of chiral symmetry. Gluon and ghost propagators in the 
presence of domains are calculated explicitly and their analytical properties 
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are discussed. The Fourier transforms of the propagators are entire functions 

and thus describe confined dynamical fields. 

Pacs: 12. 38. Aw 12.38.Lg 14.70.Dj 14.65.Bt ll.15.Tk 



I. INTRODUCTION 

Nearly every approach to the problem of the QCD ground state accepts that the vac- 
uum is characterised by strong background gluon fields and, as results of lattice calculations 
suggest, by a clustered or lumpy distribution of topological charge and action density in 
configuration space. We shall refer to such structures interchangeably as "clusters" and 
"domains". They were first observed in typical lattice gauge configurations via cooling 
or smearing algorithms which incrementally suppress quantum fluctuations by locally 
minimising or at least reducing the action density. For more recent work on cooling the 
reader is referred to ||^. On the other hand, the resulting cooled gauge fields tend to give 
rise to diminished string constant indicating a loss of confinement. An alternative way of 
analysing the underlying fluctuations of topological charge density is via the chirality of 
fermionic modes in the background of topological "lumps" , as originally undertaken by 
and re-discussed recently in Used with lattice fermions with good chiral properties such 
as overlap [^|[ or domain wall fermions the method indicates localisation of low-lying 
fermionic modes with definite chirality, the very modes responsible for the chiral conden- 
sate, for example ||10[- These results can be described in terms of the instanton liquid model 



TT[ and are regarded as evidence for instantons on the lattice. While such an interpreta- 
tion connects clusters of topological charge with chiral symmetry breaking, it says nothing 
about their relevance to confinement |jl2[. However, it might be significant that 0] have 



repeated the procedure of P with overlap fermions and no smearing and still observe strong 
localisation and definite chirality of the low lying modes. This is a piece of evidence for 
the possibility that localisation of chiral fermionic modes is due to the effective degrees of 
freedom responsible for both confinement and chiral symmetry breaking. 
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Several mechanisms of confinement have been proposed since the formulation of quan- 
tum chromodynamics. All try to realise confinement as a dual-Meissner effect, and thus rely 
on a condensation of singular gauge configurations such as monopoles and vortices |T3|,|n[. 
In particular center vortices are considered as effective degrees of freedom relevant both 



to confinement and chiral symmetry breaking WM. In general, besides the above-mentioned 



configurations characterised by topologically conserved charges there exist also topologically 



trivial domain wall singularities in gauge fixed fields [|T5l. The form in which singular fields 
occur in the gauge fixed formulation varies with the gauge choice, but their presence itself 
is most probably an intrinsically unavoidable feature of nonabelian theories, universal for a 
large variety (if not all) gauge fixing prescriptions. Consensus about this has been growing 
since the pioneering works of Gribov and Singer |[T^. This suggests that the manifestation 



of singular gauge fields is linked to the type and dimensionality of the manifold of singu- 
larities rather than to the peculiarities of their realisation within a particular gauge fixing 
prescription. 

We take as a working hypothesis that an effect of this kind can be seen in the restrictive 
influence of singular gauge fields on fluctuations in the vicinity of singularities JT^jM 



formulate a simplified model which allows one to study manifestations of this effect in vacuum 
properties and quark-gluon dynamics analytically. 

The subtleties of separating fields into regular and singular parts and the behaviour of 
regular fields at the singularities are irrelevant if one could calculate the QCD functional in- 



tegral "exactly". But these issues become crucial if one undertakes approximations [jT8[. For 
example, in gauge invariant quantities singularities due to ambiguities in gauge fixing should 
not occur. In the action such a finiteness, despite singularities in the gauge field, occurs ei- 
ther due to cancellations between derivative and commutator parts in the field strength if the 
singularity is topologically nontrivial (monopole or vortex) or due to finiteness of both terms 
separately for topologically trivial configurations (domain walls). However this cancellation 
of singularities in the action density can be destroyed by unconstrained fluctuations around 
the singular fields. Thus finiteness of the action implies specific constraints on fluctuations. 



In Ref. an example of this is considered in detail for Polyakov gauge monopoles. 

We formulate a model partition function which incorporates singularities in gauge fields 
effectively via their restrictive effect on fluctuations. We assume that singularities are present 
in general in gauge potentials, and in their vicinity one can divide an arbitrary field, A, into 
singular S and regular Q parts: 

^;:(x) = 5^(a;) + Q;:(a;). (1) 

In order that A generates finite action, it must be "close enough" to a pure gauge configu- 
ration in the vicinity of the singularity, meaning that [Q, S] = and that the field strength 
for pure gauge 5* vanishes, 

S,, = d,S, - d,S, + ig[S„ S,] = 0. (2) 

This can be realised in two ways. If Eq.(^ is satisfied via a cancellation between derivative 
and commutator parts then the singularity in 5* is topologically non-trivial and non-abelian. 
If the two parts separately vanish then the singularity is topologically trivial. The gauge 
potential S is then abelian, namely a constant unit colour vector can be associated with 
the field S. 

To be explicit, at the cost of generality, we shall take the second of these possibilities 
and further assume that singularities in vector potentials are concentrated on hypersurfaces 
dVj (j = 1, . . . , A^) in Euclidean space of volume V, in the vicinity of which gauge fields can 
be divided as above into a sum of a singular pure gauge Sj^^ and regular fluctuation part 
Q[p, with a colour vector associated with S^^\ For such fields to have finite action the 
fluctuations charged with respect to nj must obey specific conditions on dVj. The interiors 
of these regions thus constitute "domains" Vj. Demanding finiteness of the classical action 
density, one arrives at 

n^Qif) = 0, ilj = -i fe'^'^^^^l), ilj = iji fe-'''^^^ (3) 

for X G dVj, with the adjoint matrix nj = T"'n°j in the condition for gluons, and a bag-like 
boundary condition for quarks, ?7^(a;) being a unit vector normal to dVj. 
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Eqs. (1^) indicate that gauge modes neutral with respect to n"j are not restricted and 
provide for interactions between domains. In a given domain Vj the effect of fluctuations 
in the rest of the system is manifested by an external gauge field B'-^ neutral with respect 
to n'j. This motivates an approximation in which domains are treated as decoupled but, 
simultaneously, with a compensating mean field introduced in their interiors. The model 
becomes analytically tractable if we consider spherical domains with fixed radius R and 
approximate the mean field in Vj by a covariantly constant (anti-) self-dual configuration 
with the field strength 

= n^'^Bj^J, = ±Bli>, B^B^;} = B\„ 
nU) = t3 cos^^. + ^8 ^^^^ ^ ^ |(2A; + l)7r/6}to, (4) 

where the parameter B = const is the same for all domains and the constant matrix n|t" 
belongs to the Cartan subalgebra, the generators t°- being in the fundamental representation. 
Note that since the mean field represents an effect of fluctuations outside the domain there is 
no source for this field on the boundary and therefore it should be treated as strictly homo- 
geneous in all further calculations. The homogeneity itself appears as a simplifying approx- 
imation. Due to uniformity (on average) of the system outside the domain, slowly varying 
fields should be taken into account first of all, with leading contributions to this coming 
from covariantly constant fields inside and on the boundary of a domain, with a dominance 
of (anti-) self-dual fields, since they are expected to have lower action density |[T9|j20|l than 
arbitrary constant fields (see also the Ginsburg-Landau type consideration in Appendix B, 
where the appearance of a discrete set of values is also motivated). 

The model for partition function we postulate and then use for calculations describes 
a statistical system of finite density = N/V composed of ^ oo noninteracting 
spherical regions in a total Euclidean volume V —>■ oo, each of which is characterised by 
a set of internal parameters with random values: the angle uj between chromoelectric and 
chromomagnetic fields, spherical angles ipj and 6j of the chromomagnetic field, the angle 
in the colour matrix fij, chirality violating angle aj entering fermionic boundary condition 
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and the coordinate Zj of a domain. Clusters are characterised also by the fluctuation flelds 
Q^, if)^ and ip^ satisfying boundary conditions (|^), whose dynamics is driven by the QCD 
Lagrangian in the presence of the mean fleld. The propagators of fluctuation flelds for a 
given background and boundary condition can be found analytically. Thus this partition 
function offers a systematic prescription for calculation of the correlation functions, based 
on a decomposition over fluctuations and taking the mean fleld into account explicitly. Such 
a treatment of fluctuations as perturbations of a certain background fleld is sensible only if 
the essential features of the system can be seen in the lowest orders of the decomposition (at 
least semi-quantitatively). In other words one has to verify whether such basic phenomena 
as conflnement and spontaneous chiral symmetry breaking are provided by the domain- 
structured mean fleld and boundary conditions under consideration. 

In the zeroth order of the expansion we shall flnd that the gluon condensate, topological 
susceptibility x and the string constant a for colour group SU (3) take the compact form 

, 2 1^ V3 fdx . 2^3 fdx . 
j{z) = — 3 / — sm a; / — sm x , 

3z 2z J X z J X 

\ ° ; 

while the quark condensate density at the domain center, calculated in the lowest nonvan- 
ishing order over quark flelds, reads 

2 °? Jj. 2t-2(cotht-l)/2 

F{z) =e' -Z-1 + — / — 2 (cotM - 1), q = B^R^ 16, 

4 J sinh t 



where q is the absolute value of the topological charge associated with a single domain. 

To gain numerical estimates of these quantities we flxed the mean fleld strength param- 
eter B and domain radius R to flt the known value of the string tension, 

Vb = 947MeV, R'^ = 760MeV, (5) 

which leads to the values 



{as/TT){F^) = 0.081 GeV^ = 420 MeV, 

X = (197 MeV)^ {^tfj) = -(228MeV)^ (6) 

with domain charge q = 0.15 and density v^^ = 42.3fm^^. This estimation shows high 
density of clusters and strong background fields in the system, with confinement of static 
charges and spontaneously broken chiral symmetry. There is no separation of scales charac- 
terising the system, \/BR ^ 1 . The qualitative picture as well as numerical values obtained 
indicate consistency of the gross features of the model. 

These results suggest that formation of clusters, predominantly (anti-)self-dual and with 
average size 2R ^ 0.5fm, can have purely quantum origin whose explanation could require 
reference to the existence of obstructions in gauge fixing. 

It should be noted that the physical content of the above numbers can differ from other 



approaches. For instance the QCD sum rules |21| determination of the gluon condensate 
is not exactly comparable to ours, since in our case corrections of order 0{as) and higher 
contain nonperturbative information via explicit dependence of quark and gluon propagators 
on the mean field. 

Moreover, the above parameters only give a characterisation of the "bulk" properties of 
the theory and say little about confinement of dynamical colour modes and hadronisation. 
Information about these aspects is contained in the quark and gluon Green's functions, 
in particular in their propagators. It is shown that, as expected, the Dirichlet boundary 
condition removes gluon zero modes, and the propagator in this problem is well defined, 
unlike the analogous problem in the infinite volume. Both propagators have support in the 
interior of the hypersphere, where they have the usual ultraviolet singularity. Thus at short 
distances the propagators have standard perturbative form plus power corrections. The 
singularity in the configuration representation of propagators is integrable and their Fourier 
transforms exist, so that in momentum space the propagators are entire analytical functions 
due to their compact support. This we regard as a manifestation of the confinement of 
dynamical fields. 
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The paper is organised as follows. In Section 2 the boundary conditions are discussed 
and the model partition function is defined. We consider properties of the ensemble of mean 
fields in Section 3 and estimate the lowest dimension gluon condensate, the topological 
susceptibility, the string constant and the quark condensate in the lowest nonvanishing 
order in fluctuation fields. Gluon and ghost propagators are calculated in Section 4 and 
their analytical properties are discussed. In Section 5 we give an outline of the problems 
remaining to be solved and possible perspectives. The appendices contain some technical 
and illustrative material. 

II. THE PARTITION FUNCTION 

In this Section we formulate a partition function which will be used in subsequent sections 
for modelling the QCD partition function in the presence of clustered background fields. It 
should be clear from the very beginning that we will not derive the model to be considered 
from the original QCD functional integral. The mathematically accurate framework for such 
a derivation, a self-consistent mean field approximation requiring calculation of the effective 
action of QCD as a functional of the mean field and characteristic functions of the domains, 
is yet to be formulated. The best that can be done at this stage is to identify several 
ingredients of the formalism required for motivating the model within QCD, postulate the 
model partition function, and then look for signatures of justification a posteriori, by means 
of explicit calculations. 

Clustered structure of the gauge fields is introduced by the proposition that singular con- 
figurations may not be excluded ad hoc from the functional space of integration; rather the 
character of singularities should be restricted by the natural requirement that the classical 
action density for a given (in general, singular) configuration has to be finite. 

We assume that in the vicinity of singularity an arbitrary gluon field A can be divided 
as in Eqs. A^{x) = S^{x) + with Q a regular field and S the singular pure 

gauge part. 
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S,ix) = -[d^^Uix)]U-\x), U{x) = e'^f^^\ f = rt\ 



The field strength corresponding to vanishes. In this paper we will consider abelian 
singular configurations 



This can be implemented via 



N 

f = % = const, [d^,d„]fj{x) = 0, [nj,nk] = 0, 

where each of the functions fj is singular on a three-dimensional boundary dVj of the (four- 
dimensional) region Vj, while the matrices nj belong to the Cartan subalgebra of SU{3) and 
can be parametrised by 

fij = t-s cos^j + ts sin^j, < < 27r. 

The boundaries of the densely packed regions Vj necessarily intersect each other and, for 
instance, colour orientation associated with the boundary becomes ambiguous in the in- 
tersection regions. Strictly speaking, this means that the abelian singular fields should 
be accompanied by topologically nontrivial vortex-like configurations, such that the three- 
dimensional "domain wall" (a topologically trivial object) should start and finish at the 
two-dimensional singular surfaces, corresponding to a type of dislocation. A complete pic- 
ture would include the whole hierarchy of singular fields: domain walls, vortices, monopoles 
and instantons. It is hard to formulate a complete approach in a precise way. A qualitative 



discussion of this aspect of domain-like structures can be found in Ref. [Q. Even if we 
neglect the effects of "dislocations" on the boundaries dVj, a self-consistent consideration is 
still a complicated problem. However, in this case one can get some idea about features of 
the required formalism by means of an artificial example - QED in the presence of the singu- 
lar background fields, considered in Appendix A . The linearity of electrodynamics enables 
a formal definition of the free energy (effective action) as a functional of a background field 
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and characteristic functions of clusters, and thus relegates the question about formation of 
clusters in a typical gauge field configuration to a competition between energy and entropy. 
In the case of an abelian weakly interacting theory one hardly expects domain formation. 
On the contrary in non-abelian strongly interacting theory singular fields are most probably 
unavoidable, but unlike QED a straightforward formulation is a difficult task. 

First of all we should determine the appropriate boundary conditions for the fluctuation 
fields about the singular field S for finiteness of the action density. Substituting Eq. (||) into 
the QCD Lagrangian, we obtain 



:qcd = -\q%Q% + i^-m + g^ + g^ iIj 



2 
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and Q^j, is the usual field strength tensor for the fiuctuation field. We see from Eq. (|^ that 
conditions on the gluon field arise 

UjQ^ = for X G dVj, (8) 

while quark fields should satisfy the condition given in Eq. (^. Eq. (|^) means that the modes 
of the gluon field longitudinal to the colour vector are not restricted, so it is convenient 
to decompose gluon fiuctuations inside the region Vj into transverse and longitudinal parts 
with respect to n| 

Aj^^ = for X G dVj. (9) 

The separation in Eq. (|I]) into singular and regular parts imposes certain restrictions on the 
gauge transformations if the original and transformed fields Q are subject to the same bound- 
ary conditions. To determine these restrictions it is sufficient to consider the infinitesimal 
transformation 

5QI = d.u'^ - r''uj\n^B^ + AI + 5^), (10) 
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from which we conclude that gauge functions should satisfy conditions 



'1 




r — 



for X G dVj. 



(11) 



The longitudinal functions uj need not be restricted. The condition Eq.(|TT]) dictates that 
gauge fixing for the fields Q should be achieved by means of restricted gauge transformations. 



within the fc— th region with the field fluctuations in the rest of the system can be seen as 
a coupling to external gauge fields which are longitudinal to the colour direction of the 
boundary dVk- This feature motivates an approximative treatment of the partition function, 
in which clusters are treated as decoupled but, by way of compensation, a mean field is 
introduced in their interior. A self-consistent mean field approach requires calculation of the 
effective action as a functional of the mean field and characteristic functions of the domains. 
Its minima would contain information about mean field character, shape and typical domain 
size. 

Here we assume that the effective action favours formation of clusters with typical size 
R and nonzero mean field. In Appendix B it is shown that with this and an arbitrary 
constant mean field the effective action for a domain exhibits twelve degenerate discrete 
minima corresponding to (anti-)self-dual configurations and six values (for SU{3)) of the 
angle ^ associated with the Weyl group. There is also a degeneracy in the orientation of 
the chromo magnetic field. The value = is specific for an ansatz with the effective 
action polynomial in TiB'^, but the period vr/S is universal. Since the volume of the domain 
is finite the degenerate minima do not correspond to thermodynamical phases and have to 
be summed in the partition function. 

The partition function for the model is defined as 



The original conditions Eqs. (|^) show that the interaction of quark and gluon fluctuations 




(12) 
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where the thermodynamic hmit assumes V, N ^ oo with the density = N/V taken 
finite. The fields Q^^\ "ipi and tpi are subject to boundary conditions Eq.(^, in which the 
original singularities are effectively encoded. Interaction between the original domains is 
substituted by the mean field. A background gauge condition is imposed. The Faddeev- 
Popov determinant should be calculated on a restricted space of functions consistent with 
Eq. (|TT|), which can be implemented in the form of integral over ghost fields (h'^, h":) subject 
to the boundary condition 

rijhj = for X G dVj. (13) 

The integration measure dai is 

^ 2-7r 277 TT 

j dai. .. = J dai J dipi J dOi sin 9i x (14) 







2tt 5 

/ diu^ J2 SitUi-nk) I ci^,^5('^,-(2/ + l)'' 

; n 1 7 n V 



fc = 0,l 

Here yjj and 9i are the spherical angles of the chromomagnetic field, uji is the angle between 
the chromomagnetic and chromoelectric fields, S,i is the angle in the colour matrix n^, is 
the chiral angle and Zi is the centre of the domain Vi with the boundary 

{x-z,f = R\ 

The partition function Eq. (|12D describes a statistical system of density composed of 
noninteracting hyperspherical clusters, each of which is characterised by a set of internal pa- 
rameters and whose internal dynamics are represented by the fiuctuation fields. Correlation 
functions can be calculated taking the mean field into account explicitly and decomposing 
over the fiuctuations. First of all we consider vacuum characteristics of the system to zeroth 
order in this expansion. 

III. VACUUM PROPERTIES TO LOWEST ORDER IN FLUCTUATIONS 



The above prescribed perturbative treatment of fiuctuations means in particular that 
they cannot change vacuum properties of the system. Thus our immediate task is to test 
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whether the mean field itself reproduces the main nonperturbative characteristics of the pure 
gluonic vacuum. To achieve this we have to compute vacuum expectation values of a number 
of basic quantities omitting integration over fluctuation flelds. Thus we will calculate ra-point 
connected correlators of fleld strength and thereby the corresponding gluon condensates, 
string constant and topological susceptibility. 



A. Mean field correlators 

A straightforward application of Eq. ([T2|) to the vacuum expectation value of a product 
of n fleld strength tensors, each of the form 

TV 

j 

gives for the connected n-point correlation function 

(B^' (xi) . . . B'^r, (xn)) = lim y / ^ / rf(T,n(J>i . . . n(^>"Rf ) . . . B^p^ 

X ^(1 - (xi - z,f/R^) ... 0(1 - (x„ - z^f/R^) 
= BX\-.Z,,„.„^nix,,...,x^l (15) 

where the tensor t is given by the integral 

j.ai...a„ _ rlnn^^^'^^ »,0>n rO) f>U) 

''^llUl,...,^lnU„ — J "Oj't -^fllUi ■ ■ ■ -^flnfn^ 

and can be calculated explicitly using the measure, Eq. ([T^). This tensor vanishes for odd 
n. In particular, the integral over spatial directions is deflned by the generating formula 



47r J ""^^ h "^^^ j2B^u J ±j J ] ^ ^ 

where the plus and minus correspond to B'^^} being self-dual or anti-self-dual. The 
translation-invariant function 

S„(xi, ...,Xn) = lJ d'ze{l - (xi - zf/R^) ... 0(1 - (x„ - zf/R^) (17) 
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can be seen as the volume of the region of overlap of n hyperspheres of radius R and centres 
{xi, . . . , Xn), normalised to the volume of a single hypersphere v = 

S„ = 1, for = . . . = Xn- 

It is obvious from this geometrical interpretation that S„ is a continuous function and 
vanishes if the distance between any two points Xj| > 2R; correlations in the background 
field have finite range 2R. The Fourier transform of S„ is then an entire analytical function 
and thus correlations do not have particle interpretation. It should be stressed that the 
statistical ensemble of background fields is not Gaussian since all connected correlators are 
independent of each other and cannot be reduced to the two-point correlations. 

As a simplest application of the above correlators we get a gluon condensate density 
which to this approximation is 

/(^;.(^)^;.(^)) = 45'- (18) 

Note that the coupling constant is absorbed into the gauge field. 



B. Topological Charge and Susceptibility 

Another vacuum parameter which plays a significant role in the resolution of the Ua{^) 
problem is the topological susceptibility ||23|^5|. To define this we consider first the topo- 
logical charge density for the colour group SU{3), 

^2 



which in the mean field approximation takes the form 

Qi^) = ^ E ^[1 - - %)V^'] COS.;,, (19) 

where uj G {0, tt} depending on the duality of the j-th domain. We thus see that topological 
charge density is constant in each domain, and the sign of this constant is uncorrelated. For 
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a given field configuration then the topological charge is additive 

g = / d'^xQix) = q{N+ -N_), g = B^R^ /16, -Nq < Q < Nq 

J V 

where g is a 'unit' topological charge, namely the absolute value of the topological charge 
of a single domain, and i^-) is the number of domains with (anti-)self-dual field, N — 
+ A^_. With fixed total number of domains N the probability of finding the topological 
charge Q in a given configuration is given by the distribution 

""^^^ A/W 2^{N/2-Q/2q)\{N/2 + Q/2qy: 

where A/iv(Q) is the number of configurations with a given charge and Aat is the total number 
of configurations. The distribution is symmetric about Q — 0, where it has a maximum for 

even. For odd the maximum is at Q = ±g. We conclude that topological charge 
averaged over the ensemble of clusters vanishes. 

The topological susceptibility 

X = / d'x{Q(x)Q(0)) (20) 

is determined by the two-point correlator of topological charge density, which in the lowest 
approximation reads 

{Q(x)Q(y))^-^^E,(x-y), (21) 



and we get 



C. Area law for the Wilson loop 



In the same mean field approximation the Wilson loop is given by the integral 

^ r d^z- r 1 ( r - i 
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where the measure daj corresponds to an integral over the set of parameters 

{Zk, (f>k, Ok, ^k, ^k}k=l 

of the field strength 

k 

Note that path ordering in our case is not necessary since the matrices n^''^ are assumed to 
be in the Cartan subalgebra. 

Strictly speaking the contour Sl around which the path-ordered exponential is integrated 
should be a rectangle whose Euclidean-time length should be taken arbitrarily large before 
the spatial length. It is for such a contour that one has a strict interpretation of the 



behaviour of the exponent in terms of a static potential p6|j27[| . However the expectation 
that there be an area law is not dependent on the specific geometry of the contour. In view 
of the rotational properties of our approximation to the vacuum fields, it is computationally 
more convenient to consider a czrcu/ar contour in the {x3,X4) plane of radius L with centre 
at the origin. If an area law is established, as will be the case, the numerical value of the 
resulting string constant would not be precisely that corresponding to a rectangular contour. 
However due to the fact that the loop must be taken large in order to extract the potential 
the difference between a circle and a rectangle should not lead to radically different values 
of the string constant. 

To illustrate the steps in the calculation while avoiding cumbersome formulae we consider 
here the case of colour group SU{2). The details of SU{3) will be given in Appendix C, 
though the final result will be quoted below. For colour SU (2) we have 

= eVa, = ±1. 

The thermodynamic limit (V, — > oo) assumes that the subvolume 

V = V/N = 

is fixed. Calculation of the trace in colour space leads to the result 
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where we have denoted 



JSl 



(23) 



Using the properties of the measure of integration over the collective coordinates one gets 



W(L) = lim 



0) 



N 



IV V J '2 

We have exploited here the property that the integral over collective variables does not 
depend on the index j. As the contour of the Wilson loop is in the {x^, X4)-plane, the only 
nonzero components of J^j^ are 



Jm — —Ja?,\z 



f dx3dx49{l- {x- zf/R^), 



(24) 



and 



B^.J^,u{z) = 2^43 J43(^) = 2S3^43(^) = 25J43(^) COS ^, 



(25) 



where 9 is the angle between chromoelectric field E and the third coordinate axis. Now we 
can calculate the integral over the spatial orientations of the vacuum field 

' ^d^,sin^,e— -^^-^^^-^^^ 

^ 2BJa 

and the Wilson loop takes the form 



y^^^e 477 70 '^^J/^^^^'^^J'' 2BMzj) ' 



W(L) = lim 



dz 



sin 2BJ43{z) 



y Jv 2BJ43{z) 
Calculating the integral over z we obtain finally 



N 



W{L) = lim 

N^oo 



n N 



— e 



-U{L) 



U{L) 



1 



dx . \ 7rV4„3, 1 „4 

I — sm a; H -R^L + -R^ 

27iBR^ Jo X V \3 2 



7r2(l ^ cos27r5i?2) An^L rV^^^ 



v{2nBY 
Jo J(l 



+ 



;(27rS)3/2 



f 

Jo 



dx sin X 



cos- 



V Jo SL^ (2(/? — sin(/7 + s(2'0 — sin'^)) 

t-s + 1 ip t + s-1 

p — , cos — = = — , 

2Vt 2 2Vst 



(26) 
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where the thermodynamic hmit (A^, V — oo, v = V/N = has been taken. One can 

check that U{L) = when i? ^ 0, as it should. 

In the hmit of large Wilson loop L ^ R the behaviour of the first four terms in U (L) can 
be determined by inspection. For the last term a slightly more involved calculation gives 
the large L behaviour 

7r2L4 fRVL^ r{i+V-s? sin [BL"^ {2^ - sin + s(2V' - sin V'))] Stc^LR^ 
as I at - 



V Jo Ai-v^)^ BL'^ {2(p — simp + s{2iIj — sin 3v 

with corrections coming at 0{R'^). Thus only the first term in U{L), going like L^, displays 
an area dependence. The final result for the string constant for SU (2) is: 

W{L) = a = Bf{nBR^), 

For the case of SU{3), as shown in Appendix C, the function f{z) turns out to be: 

/(2;) = — 3-— / — smx / — smx . (27) 

6Z \ 2z Jo X z Jo X J 

It is positive for ^ > and has a maximum for z = l.SSvr. We choose this maximum to 
estimate the model parameters by fitting the string constant to the lattice result, 

Vb = 947MeV, i?"^ = 760MeV, (28) 

with unit charge q = 0.15, density = 42.3fm~^ and the "observable" gluonic parameters 
of the vacuum 

= 420 MeV, x = (197 MeV)^ 
(a,/7r)(F2) =0.081 GeV^ (29) 

The high density ensures area law dominance already at distances 2L «i 1.5 — 2 fm. 

The result for the gluon condensate is larger than the most recent estimate within QCD 



sum rules PSI. As already mentioned, our value is not directly comparable to sum rules 



results due to differences in the content of 0{as) corrections. What appears to be important 



is that all these quantities are nonzero, and their values can be fit to the expected numbers 
simultaneously. 

Obviously, if B goes to zero then the string constant vanishes. This underscores the role 
of the gluon condensate in the confinement of static charges. On the other hand we can 
also see that if the number of domains is fixed and the thermodynamic limit is defined as 
V, R ^ oo, N = const < oo, namely if the clusters are macroscopically large, then W{L) = 1, 
which indicates absence of a linear potential between static (infinitely heavy) charges in a 
purely homogeneous field. However this does not mean that heavy quarks {rriQ ^ B) are not 
confined if domains are macroscopically large. As is shown in Ref. [^], the nonrelativistic 



potential is quadratic in the distance between heavy quarks with the coefficient proportional 

to TUqK 

Since we have integrated over background fields exactly the role of finite range of correla- 
tion functions is hidden in the above calculation. In order to see this role explicitly one would 
need to decompose the integrand into an infinite series and integrate term by term. At this 
step all correlation functions of the background field up to infinite order would be manifest. 



The arguments of Refs. pOI about the crucial importance of a fast decay of correlators for 
confinement of static charges would be seen to apply here via this representation. 

A comment on the values of the parameters R and B appearing in our estimation is in 
order. We observe that there is no separation of the two scales characterising the vacuum. 
The average strength of vacuum fields B and the average domain size R, are comparable to 
each other \fBR 1. Neither large domains nor stochasticity of background fields are seen 
here which a posteriori justifies the mean field averaging prescription in the partition func- 
tion. This prescription corresponds to a system less ordered than, for instance, a spin glass. 



Nor does the partition function represent a heterophase mixture , since the condition for 
quasi-equilibrium is not satisfied: one may not think of these clusters as droplets of different 
thermodynamic phases as they are too small and too transient compared to the basic scale 
of interactions determined in this picture by the gluon condensate value. The mean field in 
the clusters is singled out not due to a hierarchy of scales, but due to certain specific prop- 
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erties: the (anti-)self-duality, and the abehan character if dislocations at the boundaries are 
neglected. Homogeneity of the background field appears as an approximation. 



D. Quark Condensate Density at Domain Centre 

A complete consideration of the fermionic eigen-problem for the background field and 
boundary conditions under consideration will be given in a separate work. However to com- 
plete the picture of vacuum properties in the model, we estimate here the quark condensate 
density at the domain centre. 

A complete calculation of the quark condensate in the lowest nonvanishing order over 
fluctuations requires solution of the equations 

{iP-m)S{x,y) = -5{x,y), (30) 
I fj{x)e"'^^Six, y) = -Six, y), (x- zf = R\ (31) 
S{x, y)i me-"''' = Six, y), (y - zf = R\ (32) 

where ri^{x) = {x — z)^/\x — z\, and is the covariant derivative in the fundamental 
representation, 

Substituting 

s = {ip + m)[p^no + p^o+n+i + p^o^n^i], (33) 

into Eq. (0) where 

iV± = i(l±n/|n|), 
S± = ^(l±S5/i?), 

and B = \n\B, shows that the scalar functions TY^, with ( = 0, ±1, should satisfy the 
equations: 
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(-L>2 + + 2CB)nc = 5{x, y). 



(34) 



We note that if solutions vanishing at infinity were sought, then the Green function H—i 
would be divergent in the massless limit due to the contribution of zero modes of the Dirac 
operator in the presence of the (anti-)self-dual homogeneous field. The present bag-like 
boundary conditions remove zero eigen-values from the spectrum, and the massless limit is 
regular. Due to averaging over self- and anti-self-dual configurations and all possible values 
of angle a in the partition function chiral symmetry is not broken explicitly. However, as 
we show below, a nonzero quark condensate arises in the massless limit due to an interplay 
of random distribution of the domains with self- and anti-self-dual field and the boundary 
conditions with the random value of the chirality violating angle a. 

In order to avoid cumbersome calculations and expose the role of the former zero modes 
in a transparent way we turn to the particular choice y — z — and calculate the value of 
the quark condensate at the centre of the domain. In this case the functions can depend 
only on x^, B^^^x^ and 77^ = x^/\/x^, and hence are functions of only, and the general 
solutions for scalar Green's functions take the form 



Here A(x^|/x) is the vanishing at infinity scalar propagator with mass 2B^ in the homo- 
geneous ( ant i-) self-dual field and $ is a solution to the homogeneous equation regular at 
x"^ — 0, expressed in terms of the confluent hypergeometric function. The constants can 
be used to fit the boundary condition. Terms with TYo and 7i+i are regular in the massless 
limit and cannot contribute to the trace of the quark propagator. Thus we concentrate on 
the term Using identities 



where /^^ = w? /2B + ^, and 




/^M(l + At,2,5xV2). 



7u.B^pXpP^T,+ = iB '^P^T.+, 



^lj,Bu,pXpP:^T,^ = -iB J«fPq.S_, 
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one can show that the boundary condition is satisfied if on the boundary 
which imphes that in the massless hmit C_i takes the form 

C-i = -T^2 + ir^nBR'm + 0(1), 

F{z) = e' - z -1 + — 2 (cotht-1), 

4 JO sinh t 

Moreover the singular terms cancel in pTi-i and 

hm mH_i(x,0) = ^^F(5/^V2)e-^-'/^ (35) 

and thus, 

Tr^(0,0) = ^E^(^^V2). (36) 



It should be noted that the part of the propagator responsible for the nonzero trace, Eq. (p5[), 
is proportional to the zero mode of the Dirac operator 



ipP^O^e-^'' 1^ = 0. 



Now we have to average Eg. (|36D over domain configurations taking into account the 
quark determinant. According to |^2[ the a-dependence of the quark determinant is 



exp J dxB^^B^^eil - xV-R^)| = exp{±iga} 

where q is the topological charge associated with a domain, and thus a 6'-term is generated 
effectively by the quark determinant. After averaging over a we get a finite value for the 
condensate at the center of the domain 

Numerically this is equal to {^ip) = — (228MeV)^ for B and R fixed by the string constant 
as in in Eq. (^). 
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IV. PROPAGATORS IN THE PRESENCE OF DOMAINS 



In order to study in more detail the influence of domain structure and the mean field 
on the properties of the dynamical quarks and gluons we have to find their propagators. 
They can be analytically calculated by reduction to the scalar problem, essentially that of a 
four-dimensional harmonic oscillator with total angular momentum coupled to the external 
field. The general solution is given by decomposition over hyperspherical harmonics. In the 
following section we present the exact solution for the scalar propagator, though with most 
derivations relegated to Appendices D and E. With the scalar result we derive propagators 
for ghost and gluon fluctuations in an external (anti-)self-dual field with Dirichlet boundary 
conditions imposed on the fluctuations on a hyperspherical surface. 

A. Scalar Propagator 

The problem to be solved is given by the scalar Green's function equation 
[{d, - iB,f - Ae]G{x, x'\ix) = -5^'\x - x'), 

B/u = --Bi^i,x„, B^^Bfj^p = B'^S^p, (37) 
with the homogeneous Dirichlet boundary condition 

G{X,X'\I1)^2^R2 = G{X,X'\IJ)^'2^R2 = 0, 

where R is the radius of a hypersphere centred at the origin and ji — M?" jlB. 
We present first the solution to the corresponding eigenvalue problem. 

One may choose B^^ such that 

S34 = £/, B\2 — -B, — B ^ -£/ < B. 
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A representation of the eigenfunctions in terms of a complete orthonormalised set of eigen- 
functions of the four-dimensional Laplace operator is achieved in the following hyperspherical 
coordinate system (see for example, p3| ) 



Xi = r sm 1] cos (p 
X2 = r sin rj sin 
X3 = r cos T] cos X 

X4 = r cos T] sin x- (38) 



The angular eigenfunctions are 

C,^,mM<P,x) = (-l)l™^+"^l(27r)-ier-'"^'™^+"^(r/), 
X exp z[(mi — m2)x + ("^1 + "^2)0] 



0fc"'"'''"'(^) = ^J2{k + l){k - ry.{k - s)\r\s\ 

^ (-l)'-"cos'=-'^-^+2"r7sin^+"-2'"?7 
(/c — r — s + n)!n!(r — n)!(s — n)! 
(r, s = 0, 

where fc, mi, m2 are respectively the orbital angular momentum and the two azimuthal quan- 
tum numbers, relevant for a four- dimensional hyperspherical symmetry. That the Ckmim2 
are eigenfunctions with the said eigenvalues is proven in Appendix D. 

The eigenfunctions for the complete problem are a product of radial and angular parts. 

The radial equation has a solution expressed in terms of the confluent hypergeometric func- 



tion, in the notation of \M 



fir) = {Br'/2f/'e~^^'"M (^- + 1- m,,, + k + 2; Br' /'^ , (39) 

where the function regular at r = is chosen for normalisability. Here m2,i should be put 
equal to m2 for the self-dual field, and nii for the anti-self-dual field. It is convenient to 
denote 
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k 

n2,i = - + l-m2,i. (40) 

Another independent solution (not normalisable in our problem) , which is regular at infinity 
and singular at the origin, would be obtained by replacing the function M with the function 
U. Imposition of the Dirichlet condition aXr = R forces eigenvalues A to take discrete values 
defined by the zeroes of M(a, 6, z) as a function of a at fixed h and z. The eigenvalues A are 
strictly positive. As we will see below, the case of M.^ = —2B will be met in the problem 
for gluon propagator. In this case the lowest eigenvalue Aq is defined by {k = mi = m2 = 0) 



m(-Ao,2;S/?V2) =0, 



and is a positive function of BR^, as can be checked. 

The propagator can be found by the standard method of decomposition over hyperspher- 
ical harmonics, see for instance |^ , exploiting the above hyperspherical representation. The 



derivation is given in Appendix E but it consists of essentially two steps. First we exploit 
completeness of the angular eigenfunctions in order to reduce the four dimensional problem 
to a one-dimensional Sturm-Liouville problem with known exact solution for the radial de- 
pendence. Second we add a solution to the homogeneous equation with coefficient selected 
to implement the above boundary condition. 

Using the solution to the eigenvalue equation with A = 0, one gets independent solutions 
for the homogeneous equation. In the notation of , the two solutions (respectively regular 
at the origin and at infinity) are 

Riir\k, n2,i, fi) = r'=e-^'''/^M(n2,i +fi,k + 2; Br^ /2) 

R2ir\k, 712,1, 1^) = r'=e-^'''/^f/(n2,i + /i. A; + 2; 5rV2). (41) 
The Sturm-Liouville equation is satisfied by 

X,.,,(r,rV) = B^^^j^^Rkn,Ar,r'\i^) (42) 

with 
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Ri{r\k,n2,i, fJ,)R2{r'\k,n2,i, fx), r < r' 
Ri{r'\k,n2,i, iJ,)R2{r\k,n2,i, IJ,), r > r' 
In terms of these quantities, the Green's function is: 



(43) 



k,rn\;rn2 



V{k + 2) 



(44) 



X 



' . /, ^ U{n2^i+ii,k + 2-BR^/2) ^ ^ ' 
Rkn2,Ar,f m - YT/ \ r~r^T7~51yri^^n'^\^^'^2,i,IJ')Ri{r \k,n2,i,li) 



M(n2,i + /i,fc + 2;5i?72)' 
The first term inside square brackets guarantees G to be a Green's function through the 
solution to the Sturm- Liouville equation. The coefficient in the second term is determined 
by the boundary condition. The only singularity in G is the usual ultraviolet one at a; = x'. 
For more details we refer the reader to Appendix D. 



B. Ghost and Gluon Propagators 



We work in the background gauge 



D'^'Al = 0, 



and use conventions for the adjoint representation of colour SU (3) 

n^T^cos^k + T^sinCk, Cfe = (2A; + l)7r/6, A; = 0, 1, . . . , 5. 
B^u = nB^,, B = V^B. 

The colour matrix n can be diagonalised by the unitary transformation 

n = C/nC/t - diag (Ci, -Ci, 0, C2, -C2, Cs, -Cs, 0) 

Ci = sine, C2 = (sine + V^cosO/2, Ca = (- sin^ + V3cosO/2, 



where U = diag(lV, 0, W, W, 0) with 



W 



V2 



(45) 



(46) 
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For the values of the angle as in Eg. (|45|) the Q take values from the set (±l,±l/2). 
Namely, for k = 0,1,2,3,4,5 respectively, the three are (1/2,1,1/2), (1,1/2,-1/2), 
(1/2,-1/2,-1), (-1,-1/2,1/2), (—1/2,1/2,1). The diagonalized covariant derivative takes 
the form 

= UD^U^ = d^- inB^. (47) 

In the Feynman gauge, diagonalized equations for the ghost and gluon propagators take the 
form 

~D'^G{x,x') = 6{x-x'), 

(-D''5^, + 2inB^,) G^p{x, x') = S>''5{x - x'), (48) 

with the diagonalized boundary conditions 

nGup{x,x) = 0, nG{x,x) = 0, for x^ = or x'^ = R^, 

where n and the propagators are diagonal matrices in colour indices. As evident above, 
the matrix n has two zero eigenvalues and the corresponding gluon components are not 
restricted by above boundary condition so the equations for these modes are simply the free 
ones. These modes are thus not confined in the model under consideration. 

The scalar equation for the ghost propagator has been solved in the previous section, 
where one should simply replace B nB and put = 0. The equation for the gluon 
propagator and the boundary condition can be further diagonalised with respect to Lorentz 
indices and is thus reduced to four scalar equations, each of which has a well defined solution 
since, as discussed above, dangerous zero modes, so called chromons [^,^], do not satisfy 



Dirichlet conditions and do not contribute to the Green's function. The original propagators 
are restored by the inverse transformations. 

In the case of the gluon propagator one can avoid the second diagonalization by looking 
for a solution of the form 

G,p{x, x') = [d^6,p + 2tnB,p) A{x, x'). (49) 
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Substitution of Eq. (^91) into the original Green's function equation gives 

-{D^Y + An'B^] A(x, x') = 5{x - x'), 
with the solution given formally by 



1(1 1 



6{x~x'). (50) 



A\n\B y-D"^ + 2\n\B) -D^-2\n\B^ 

The two terms above are nothing but scalar propagators with "mass term" Ai^ = ±2\n\B. 
Substituting Eq. (|^) into Eq. (^9|) and using notation of the previous subsection for the 
scalar propagator with B —>■ \n\B and n = ±1 one gets after simple manipulations 

G,,{x,x') = -V [Gix,x'\l) + G{x,x'\ - 1)] + — ^ [G{x,x'\l) - G{x,x'\ - 1)] , (51) 
Z Z\n\JD 

where only nonzero elements of the diagonal matrix n are involved. With this representa- 
tion it is clear that the boundary condition for the gluon propagator is satisfied if the scalar 
Green's functions G{x,x'\ ± 1) are subject to the homogeneous Dirichlet condition indepen- 
dently of each other. An explicit form is obtained via Eq. (^^ by substitution B \n\B 
and /i = ±1. This can be done straighforwardly for the terms with k — m2,i > in the 
expansion over hyperspherical harmonics, as is obvious from the integral representations of 
the confluent hypergeometric functions [0 (6 > a > 0) 

M(a, &, z) = ^ /' dte^H'^~\l - tf-^-\ 

1 (a)l (o — a) Jo 

1 

U{a, b,z) = -— dte-'Y-\l + tf-''-\ (52) 
1 \a) Jo 

with a = 1 + fi + k/2 — m2,i and 6 = A; + 2, in our case. Special comment is required for the 
terms with a = k/2 — m2,i = in the decomposition of G{x,x'\ — 1). 



Using the representations Eqs.(52) it immediately follows that 



lim M(a, b, z) = 1 + 0(a), lim U (a, b,z) = l + 0{a). 

a^O a^O 
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With this and Eqs. (^lD-(^4D one can be convinced that the singularity in the gamma- 
function in Eq. (^) at a = ^2,1 — 1 = kjl — vtii^y = is cancelled by the contribution coming 
from the expression in the square brackets. Thus we conclude that the gluon propagator 
exists. This confirms the absence of the zero modes under the imposed Dirichlet boundary 
conditions. 



C. A comment on analytical properties 



The scalar propagator Eg . (^41) which determines the analytic properties of the off- 
diagonal components of the ghost and gluon propagators have compact support in hyper- 
spherical region of radius R in Euclidean space-time with the usual ultraviolet integrable 
singularity at x' = x. Thus the Fourier transform of the propagator averaged over domain 
position, given by the integral 

G{p^) = [ d^xe'P^'Gix), 
G{x — y) = dzG{x — z,y — z), (53) 

J 

leads to a G'(p^) which is an entire analytical function in the complex p"^ plane. Entire 
propagators are typical for nonlocal field theories and has been interpreted as confinement 
of dynamical charged fields |1I9|J36| , |37|| . Thus the presence of domains maintains confinement 



of off-diagonal gluons and ghosts. 

An instructive example is given by a toy calculation which illustrates the qualitative 
behaviour of the Fourier transform of propagators with compact support in a finite region 
of R'^. We calculate the Fourier transform of the function 



The calculation proceeds via the following steps (where p 

1 . , ^ . 2 /•! 



D{p) = - f dWl - I drre'P^^ = - f dtVl - f drrcos{prt] 

TT J-1 Jo 71 JO JO 
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dxJi{x) = J2k+2{Rp) 



1 [Rp 



P\=o 



1 

- / drJiipr) 
p Jo 

1 - MRp) 

9 ' 



where we have used the identity 

oo 

Jo{z) + 2Y,J2k{z) = 1. 

fc=l 

This propagator is an entire function with the properties 

loiRp) - 1 



D{0) = R'/A, D{ip) 



hm D{p) = — 

p2_»oo p^ 

hm D{ip) 



TrRp 



cos{Rp — 7r/4) 



(54) 



(55) 



which indicate the standard behaviour for asymptotically large Euclidean momenta, and 
an exponential rising in the physical region (large energy). We intend to consider elsewhere 
detailed analytic properties of the Green's functions in the present approach (including that 
for fermions). 



V. CONCLUSIONS AND OPEN PROBLEMS 

The idea of domains in the vacuum is not a new one and various hints and attempts at 
implementation of such an idea can be found |T^j38| -^. These approaches assume explicitly 
or implicitly that the boundaries of domains are populated by the (chromo)electric and/or 
(chromo) magnetic "charges and/or currents" which produce nonzero field strength inside 
domains. Thus the source for the mean field inside is assumed to be present on the boundary. 
Specific configurations suitable in principle for a description of such domains are known (see 
for instance [p^,p]-P^). In this picture domains are assumed to be stable and in this sense 



are somewhat similar to the usual domains in ferromagnets. 

The model presented in this work differs cardinally from this picture. The central idea 
that enables us to introduce and consider domains is the observation made in Refs. |1T7|JT8[| 
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that the presence of singular pure gauge background fields imposes specific conditions on 
quark, ghost and gluon fiuctuations. The boundaries correspond to the locations of singu- 
larities in the pure gauge vector potentials which by themselves do not generate any field 
strength. Such boundaries make their presence felt only via their impact on quantum fiuctu- 
ations. The mean field inside domains appears as a collective effect of quantum fluctuations, 
which themselves remain subject to certain boundary conditions. The domains are not stable 
in this picture, but describe a speciflc class of field fiuctuations in the system. Within this 
model all the fundamental features of the QCD vacuum - gluon condensation, topological 
susceptibility, confinement of static and dynamical charges and a non-zero quark condensate 
- emerge in a transparent and simple way. 

So far we have discussed this mechanism in a purely qualitative manner and the rela- 
tionship of the model with real QCD has to be clarified. It should be recalled that our 
motivation skipped over two points, both requiring more formal justification. In the first 
step we prescribed a particular way of dealing with singular pure gauges and thus the QCD 
functional integral incorporated densely packed interacting domains. In the second step 
we replaced this integral by a model partition function describing decoupled hyperspherical 
domains or clusters. The role of the mean field inside domains is to compensate effectively 
for the decoupling. The problem of verification of both steps remains open. In particular, 
a possible relationship between this kind of domain formation via singular pure gauges and 
the Gribov problem has yet to be understood. 

Concerning phenomenological applications, a complete solution to the fermionic eigen- 
value problem would immediately enable clarification of the connection between the picture 



of spontaneous chiral symmetry breaking in this model and the Banks-Casher relation [jTO 

With quark and gluon propagators in the mean field detailed applications to hadron 
physics are accessible. The meson spectrum, for example, can be computed via a bosonisa- 



tion procedure as applied in p^Ji5|] or via Bethe-Salpeter equations. Entire quark and gluon 



propagators are expected to give rise to the Regge character of the spectrum of relativistic 
bound states [p|,^. 
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In this context the Ua{^) problem can be also addressed. Preliminary estimations show 
that due to nonzero topological susceptibility the pseudoscalar correlators in the isovector 
and isoscalar channels are different in the massless limit and strong splitting between the 
masses of the rj' and tt mesons is expected, Alternately, the anomalous Ward identity of 
could be studied order by order in the decomposition over fluctuation fields. The fact 



that for the pure glue theory a reasonable value for the topological susceptibility is obtained 
simultaneous with a non-zero quark condensate is encouraging in this respect. 

It would be tempting to look for the present picture in lattice simulations. However, 
domains of constant field can only be taken seriously in a statistical sense, so one should 
compare results not configuration by configuration (say, after moderate cooling) but for 
correlators and condensates calculated within the model and on the lattice, where a full 
statistical ensemble has been taken into account. 

Returning to gluonic fluctuations, the picture of dynamical confinement remains incom- 
plete. Diagonal or "neutral" gluons remain freely propagating modes in the pre-mean field 
framework. Intuitively it is clear that this problem is ultimately related to the topologi- 
cal triviality of the class of singular field we have considered here. Incorporating a wider 
hierarchy of singular fields can resolve this problem. 

Finally, a set of open problems relate to the general properties of quantum field theory 
with domain-like structures and Dirichlet boundary conditions on fluctuation fields. Entire 
propagators, which appear as a result, indicate that the theory is nonlocal. The ramifications 



of non- locality need to be investigated, particularly in light of recent work by Efimov [47 



Although the choice of boundary conditions is not expected to generally influence short- 
distance singularities in Green's functions, the explicit structure of ultraviolet divergencies 
and the question of renormalisability of a quantum field model with Dirichlet boundary 
conditions imposed on fields in regions of space should be investigated explicitly. 
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APPENDIX A: QUANTUM ELECTRODYNAMICS 

The purpose of the following is to illustrate how the effective action as a functional of the 
mean field and characteristic functions can be defined formally in the abelian case. Consider 
QED 

L^-^Q^^Qf,^ + iP{i^-m-e(Q)ip-eip^ip (Al) 
in the presence of an external pure gauge singular field of the form 

N 
3 

where the functions fj have topologically trivial singularities on hypersurfaces dVj and are 
assumed to be not Fourier transformable. Gauge transformations which would remove such 
a pure gauge field are then not defined. Here Q^,, — d^Qi, — dj^Q^ is the field strength 
for the photon fields. Thus S appears only in the interaction term coupling to the fermion 
field. The fluctuation flelds Q, and il) are assumed to be regular differentiable functions 
everywhere in Euclidean space. It should be stressed that unlike non-abelian theory there 
is no internal necessity for considering singular fields in electrodynamics, and the example 
below is artificial in this sense. 

The field S in the vicinity of the j— th singular surface can be represented as 

S,-v'Md.)ax). (A2) 
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where is a unit vector normal to the surface dVj. Finiteness of the action density thus 
requires that 

i^ix) f{x)i;{x) = 0, X E dVj. (A3) 
This condition is satisfied if, for x on the boundary 

tlj = -I fe^'^'^^tlj, ^ = ^i fe-''^^^^, (A4) 
which is the well-known bag-like boundary condition Note that we are working in 



Euclidean space-time and the fields ip and ip are independent variables. The angle aj is 
arbitrary and need not be the same for different j. It should be stressed that the boundary 
condition violates chiral symmetry. No conditions on the photon fluctuation field Q arise 
since it is decoupled from S. Now we can write down the functional integral straightforwardly 

Z[S] = lvQ5{dQ)J^ PV^^^/>e-5[^+^'^''^], (A5) 

where now the space J-'s contains only those fields which satisfy the boundary conditions 
Eq.(^). We stress that the field S in Eq. ( [A5[ ) is considered as a fixed background field. 
Gauge fixing for the field Q can be achieved by regular gauge transformations. We see 
from Eq.( |A5D that due to the presence of the singular field S, the integral over fermionic 
fluctuations is separated into integrations over fields inside subregions Vj bounded by the 
surfaces dVj where the background field is singular, and these fluctuations are subject to 
the boundary conditions (|A4|) . 

Now we define a procedure for averaging over singular configurations. This is done by 
identifying the set of different singular configurations with a set of characteristic functions 
dividing Euclidean space into subregions, whose boundaries coincide with the singular sur- 
faces of a given singular field 



N 

{S} « {Xi, • • • , Xiv} , E y d'xx]{x) = V, 



where the requirement of conservation of the total volume is imposed. Integration over S is 
defined as an averaging over an ensemble of characteristic functions and angles aj coming 
through the fermionic boundary conditions, 
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VSZ[S] « n/ Vxjda,6 (l -V^'J^J d'xxkix)) 

j \ k / 

x2[xi,---,XN\ai,...,aN], 

N _ 

Z[x\oi] = liJVQSidQ) J Vi^jV^j 

1 



(A6) 



X exp 



/ 



(A7) 



In this representation translation invariance as well as chiral symmetry (for m = 0) are 
restored because of the averaging over all aj and characteristic functions Xj- 

Let us integrate out the photon field and factorise the part of the fermionic integral 
corresponding to A;— th region. We obtain 

_ N 

Zk[x\a\^ j ViPk'D^kU J VijjVijj 



X exp < 


Jvk 






X exp < 


k 















where Jj{x) denotes the electromagnetic current and D{x — y) is the standard photon prop- 
agator. Inserting 'unity' represented as 



1= n It^BH B^-eY^f d'yD(x-y)Jj{y) 
we arrive at the representation 
Z 



[X\a\ = j VB^e^v{-SAB^\x.(A} (A8) 
j ViJk'D^k exp I d^x^Pk{i ^ - m + e 



2 d'^xd^yJk{x)D{x - y)Jk{y) j- , 
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N _ 



X exp Ij d xipj {i ^ — m)ipj 



o2 



+ ^ E / d'x [ d'yJ^{x)D{x-y)Jy{y)\. (A9) 

In this representation the partition function Z[x\oi\ is defined by the fluctuations of the 
fermion field in an arbitrarily chosen subregion Vk in the presence of the electromagnetic 
field -B^, the dynamics of which are governed by the effective action Sgfj. 

As is seen from Eq. ( |A9|) , this effective action has appeared as an integral (or collective ef- 
fect) of field fluctuations in the rest of infinite system, outside the k—ih domain. This action 
functionally depends on a division provided by a particular set of characteristic functions. 

Physically different situations would arise depending on the properties of the effective 
action. This becomes obvious if we write down the partition function averaged over an 
ensemble of the characteristic functions and boundary conditions (angles aj), 

N , / N 



2 = n / T^X,da,5 U ~V-^Y.j d^^X^{x)\ j VB^ exp {-5eff[fi' 

j \ i / 

+y d\d%Ju{x)D{x - y)Jk{y)^ . (AlO) 

Two qualitatively different pictures are possible. If the functional S^a has an absolute 
minimum at -B^ = and infinitely small volumes of all regions Vj {j ^ k), then we recover 
the standard QED partition function in the infinite volume. If, however, the minimum is 
at nonzero mean field and some nonvanishing averaged size of subregions is supportable, we 
depart from standard electrodynamics. In principle the effective action can be calculated 
(at least within perturbation theory). As mentioned in the introduction, there is no reason 
to expect that the second scenario is realised in electrodynamics. 



36 



APPENDIX B: EFFECTIVE POTENTIAL FOR THE CONSTANT FIELD 



Consider a covariantly constant abelian field with the field strength parametrised as 



= n"r« = cos e + T« sine, 

Ei — B^i, Hi — -eijkBjk, + = 2B^, 
EH= |E||H|cosw, {EUf = U^{2B^ -H^)cos^u, 



and let the gauge invariant effective potential be given by the series 



Ues{B,uJ,0 = Y^AkTrB^'', {B''),, = n^^^, . . . B, 



(Bl) 



with Ak constants. One can show that if Ucs is bounded from below and has a nontriv- 
ial minimum as a function of parameter B then there is a set of twelve discrete minima 
corresponding to an (anti-)self-dual fields and six values of the angle ^. 

The odd powers of n and B do not appear in the potential, since this would mean 
violation of Weyl symmetry and parity respectively. The Weyl group is a discrete subgroup 
of global SU (3) and in this case can be seen as the group of permutations of the eigenvalues 
of the matrix n. Such permutations can be arranged by a shift of the angle parametrising 
the abelian field configuration, ^ ^ ^ + 7rn/3. In other words, the effective potential is 
periodic in ^ with a period 7r/3, the angle of 5'C/(3). It is also periodic in uj with the period 
TT due to invariance under parity. This can be checked using formulae 



Trn^ - 3, Trn^ = (9/4), Trn^ = (3/16) [10 + cos(60], 

TrS^ = -2(E2 + H^) = -4^2, 

TrS^ = 2[(E' + H')' - 2(EH)2] 

TrS^ = -2(E2 + h2)[(E2 + H2)2 - 3(EH)2] 

B^--CEKf 
4 



5^ - -(EH)' 

2^ ' . 
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A nontrivial dependence on ^ appears for A; > 3. Higher terms depend on ^ via functions 
cos 6/^ {I > 1). Taking into acount the first three terms in the decomposition and calculating 

the traces as above one gets 



The coefficients C2 and C3 are assumed to be positive to provide for the boundedness of 
the potential from below, and A is a scale. The sign of the constant Ci is of particular 
importance. If Ci is negative the minimum is trivial B — 0. For Ci positive potential has a 
minimum at nonzero B. Using the identity 



it is easy to check that there are degenerate absolute minima corresponding to field config- 
urations with the parameters 



These twelve discrete degenerate minima correspond to self-dual and anti-self-dual field 
configurations and six values of angle ^, and there is a continuous degeneracy relating to 
orientations of the chromomagnetic field H. With the symplest polynomial form for Ues as 
above we have ^0 = ^r/G. This value depends on the form of effective potential, however the 
period 7r/3 related to the Weyl symmetry is universal 

APPENDIX C: THE WILSON LOOP FOR SU{3) 
For SU(3) the eigenvalues of the colour matrix 



C/eff = -C,B^ + § - ^(EH)^] + ^,B\10 + cos 60 [b' - ^(EH)^' 



(EH)2 = h2(252-H2)cos2cu, 



= B^ ujn = 7rn (n = 0, 1), = {2k + l)n/6 (A; = 0, 1, . . . , 5), 




fij — cos Q + sin 



take values from the set 
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Ill 



1 



for all different values of the vacuum angle 

G((2/ + l)7r/3)z=o,...,5. 



This leads to the following result for the trace averaged over the vacuum angle 



W(L) = lim 



6 ^ 



TV 



Then the Wilson loop takes the form 



WiL) = lim 



1 U(L) 



N 



-U{L) 



U(L) = 3 — — / — smx — — / — smx 

^ ^ 3v \ 2nBR^ Jo x nBR^ Jo x j 

f V3 2 

7r2(l - cos 2t:BR^/^ 2n^{l - cos nBR^/Vs 



;(27r5/v^)2 



v{nB/Vsy 



+ 



+ 



(27r5/v^)3/2 



dx sin X 



(7r5/v^)3/2 Jo 



dx sin X 



ds I dt\ V3 



V Jo 
sin 



4^3 



(i-v^)2 \ (2v9 — sin (y9 + s(2'?/' — sin ■?/')) 

BL^ {2ip - sin + s{2ip - sin ip)) /2V3] \ 
BL"^ {2(p - sirup + s(2^ - sin^)) ) ' 



This leads to the string constant as given through the function as in Eq.i 



APPENDIX D: SCALAR FIELD EIGENVALUE PROBLEM 



We introduce the 0(4) generators 
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Mi = i{x^di - Xid^), 

respectively for spatial rotations and Euclidean "boosts". These satisfy the usual commu- 
tation relations. 

Now in the scalar field eigenvalue problem we encounter the structure 
Given that 

B^d^ — '^B^i/Xiid^ 

= ^{EM^-BL^) 

we see that it is better to go over to the 0(3) x 0(3) generators, 

Ki = i(L + M) 
K2 = ^(L-M). 

So with the field self-dual/antiself-dual, E — ±B we obtain, 



iBK2„ E^B 
iBKu, E^-B 

Then the four dimensional Laplace operator can be written 



The usual considerations show then that the eigenvalues of the complete set of mutually 
commuting operators = K|, Ki^ and are 

^(^ + 1)(^ = 0,1,...), m„m, = |^-l,...,-^. (Dl) 

Eigenfunctions corresponding to these values can be found such that 
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2z^ kmim2 



■m2Ci 



km\m2 ■) 



(D2) 



In the hyperspherical coordinates 



4K.I takes the form, 



4K? 



Xi — r sin rj cos (f) 
X2 — r sin rj sin (p 
xz — r cos T] cos X 
Xi — r cos 77 sin x 



Q2 Q2 Q2 Q 

■— - cosec'ri— - sec'r)— + (tan ,7 - cot 77)-. 



Thus 



4:KlCkmim2{'n,(l>,X) = expi[(mi - m2)x+ (mi + m2)0] 



X 



drf 



+ (mi + m2) cosec 77 + (mi — m2) sec 77 



+(tan77 — cotr?) — 

Of] 



0mi — m2 ,mi +m2 



Up to normahsation, the angular eigenfunctions © can be better written in terms of the 
hypergeometric function, 

Qmi-m2,mi+m2(^) OC COS"^^""^^ (77) sin'=-"^i+"^^ (77) 

k k 

X2Fi(-- + mi, -- - m2; mi - m2 + 1; - cot^ 77). 

For compactness of notation we denote 

h{7]) = cos"^i-'"2(77) sm''-^'+"''{r]), 
k k 

u(r]) = 2Fi(-- + mi, -- - m2;mi - m2 + 1; -cot^77). 



Then, after some tedious calculation, one gets 
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4KlCkmim2{vA,x) = expii[(mi - m2)x + ("^i + m2)0]/i(?7) — — 

arj 



+ cot ( ^ '2 - (2('^ - ^1 + "^2) + 1) 



c/m / 2(mi — + 1 
drj \ cot^ ?7 

+u ({2mi — k){2m2 + k) cot^ rj + 2k{mi — 1712 + 1) + Amim2 



Rewriting in terms of the variable z = — cot^ rj one eventually brings this to the form 

4:KjCkmim2{v,^,x) = expii[(mi -m2)x+ {mi + m2)0] x 4/1(77) (1 - ^) 



X 



~ ^^~d^ ~^ 'dz ^^^^ - 'm2 + I) - {I - k + mi - m2)z) 



■2 ''2 ^ 2'2 



(D3) 



But u satisfies the hypergeometric equation. Thus the first three set of terms in the square 
bracket of Eq.( |D3| ) vanish, leaving 



namely, Eq. ( P2| ), which completes the proof. 

Putting all this together we arrive at the following representation for the square of the 
covariant derivative 

{d^ - iB^f = ^drir^dr) - ^Kj + 2BK2, - -^B^r\ E = B 
[d, - iB,f = y^drir^dr) - ^Kl + 2BKi, - -^Bh\ E = -B. 

Using the eigenfunctions Ckmim2 we can reduce the original eigenvalue problem to that 
corresponding to a radial operator 



- [y,dr{r^dr) - + 2m2,iS - ^ - MV{r) = A/(r) 



where m2,i is defined in the main text. 

So the complete eigenfunctions will be a product of radial and angular parts, 

l/j(x) = f{r)Ckmim2{V,(l>,X)- 
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The radial function can be solved, as described in the main body of the paper, leading to 
the solution given in Eg. (p9|) . 

A comment at this point on the half-integer values of the azimuthal quantum numbers. 
The angular eigenfunctions depend only on the sum and differences of ms, which will be 
whole integers. Also the combination 

iXi th hj r\j th th rv 

n,,,^--m,,, + l = - + - + !,- + -,...,--- + ! 
= k + l,k,...,l. 

This combination, and thus the eigenvalue spectrum, will always involve integral values. 

APPENDIX E: DERIVATION OF SCALAR FIELD PROPAGATOR 

We represent the delta function in the hyperspherical coordinates, 

^(4). _ ^ 6ir-r')5iv-v')5i<l>-<l)')5ix-x') ^ 

sin 1] cos r] ' 

where the primed variables correspond to the hyperspherical coordinates of x'. We use the 
ansatz, 

G{x,x'\fi)= Vkmrm2{v\(l>',X)Ckmrm2iV,'P,X)Xkmrm2ir,r'\fi). 
fc, mi, 7712 

Inserting this into the original Green's function equation, exploiting the delta-function rep- 
resentation and the fact that the functions C are eigenfunctions of the covariant derivative 
squared operator, yields 

k,mi,m2 



X 



3 " J.2 ■ - ----- 4 

(5(r - r')6{7] - - ^')6{x - x') 

sin 1] cos rj 

We can now separate the angular from the radial dependence in two equations, 



43 



k,m\,m2 



k(k + 2) „ 5V .,2' 



sm 7] cos ?7 

(r, r'l/x) 



(El) 



(5(r — r') 



(E2) 



Equation ( |E1| ) we recognise as the completeness relation for the angular eigenf unctions, thus 
V = C. Using Eq.(^), we read off that X does not depend on both quantum numbers 
1712, nil but on one of them, depending on whether the field was self-dual or anti-self-dual. 
In fact we shall indicate this dependence on m2,i via the quantum number ^2,1 so that 
X = Xfc„2_i(r, r'lfi). 

We can solve the radial problem by solving for the radial Green's function in the infinite 
volume and then adding a solution to the homogeneous equation with an arbitrary coefficient. 
The coefficient is fixed by imposing the finite boundary condition at r = i?. 

Using our solution to the eigenvalue equation we can easily extract homogeneous solu- 
tions, namely to the equation: 



k(k + 2) ^ B^r^ .,2' 
^ ^ +2m2iB ; 



R{r) = 0. 



In the notation of the two solutions (respectively regular at infinity and the origin) are 
Ri{r% n2,i, /i) = r'=e-^'''/^M(n2,i + /i. A; + 2; Br" jl) 

R2{r% n2,i, /i) = r'=e-^'-'/V(n2,i + /x. A; + 2; Br^ jl) (E3) 

with n2,i defined in Eq. (^) . 

We next solve for the Green's function in the R = 00 case by recasting the problem in 
the form of a Sturm-Liouville equation, 

[drp{r)dr + g(r)] X(r, r') = —6{r — r'), 



so that 



p(r) = 
q{r) = - 



B^r^ k{k + 2) 



+ 



2Bm2 1 + M"" 



AA 



The Sturm-Liouville equation is known to have solution, 

1 



X{r, r') 



p{r')w{r') 

with w{r) the Wronskian of the homogeneous solutions, 



-i?(r,r') 



w — R1R2 — R2R1 



and 



R{r, r') = < 



Ri{r)R2{r'), r < r' 
Ri{r')R2{r), r > r' 
The Wronskian for the two solutions is evaluated to be 

^{k + 2)^/2B 



w{r) 



so that 



r(n2,i+/i)(5rV2)3/2 
Ar{k + 2) 



p{r)w{r) — 



BT{n2,i+ii) 

We now construct the full Green's function and impose the boundary condition at finite 



R: 

G{X,X'\^)= Ckmim2{v'A',X')CknHm2{V,(l>,X)[Xkn2,i{r,r'\^) 

+Akn2,i{R\B, ii)Ri{r\k, ^2,1, ii)Ri{r'\k, ^2,1, //) 
Imposing G{R, r'|/x) = with r' < r = R we have 



Akn2,AR\t^) 



Ri{R\k, n2,i, n)Ri{r'\k, ^2,1,/^) 
w{R)p{R)Ri{R\k, n2,i, /i)i?i(r'|/c, n2,i, /x) 

722,1, lj)R2{R\k, 722,1, /^) 

w{R)p{R)Ri{R\k, ri2,i, /i)i?i(r'|/c, rz2,i, /x) 

i?2(-R|A;, 722,1, A«) 



Thus finally. 



w{R)p{R)Ri{R\k,n2,upL) 
r(722,i + /i) f/(n2,i + /X, A; + 2; BRy2) 



(E4) 



4r(fc + 2) M(7i2,i + /X, + 2; 5i?V2) ■ 
It is straightforward to see that imposing G{r, i?|/x) = 0, r < R will give the same result for 
A. 
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